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The bosonization technique to describe correlated electrons in a one-dimensional quantum wire containing impurities is 
extended to include an applied voltage source. The external reservoirs are shown to lead to a boundary condition for the boson 
phase fields. We use the formalism to investigate the channel conductance, electroneutrality, and charging effects. 
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PACS numbers: 72.f0.-d, 73.40. Gk 



The puzzling physical properties of one-dimensional 
correlated fermions at low temperatures can most conve- 
niently be described within the bosonization technique. 
This method allows for an exact treatment of 

Coulomb interactions that hamper most other theoreti- 
cal approaches. In this Letter, we describe how a voltage 
bias can properly be incorporated in terms of a boundary 
condition for the bosonized phase fields. The formalism 
is similar in spirit to Landauer's approach developed for 
noninteracting electrons Q . The capacity of our concept 
is demonstrated for a single impurity embedded into the 
correlated one-dimensional (ID) electron liquid, a prob- 
lem that has attracted considerable theoretical interest 
[5,[|, and is beginning to find experimental realizations 

13- 

In order to describe Coulomb interactions adequately, 
one has to specify the setup under consideration. If one 
deals with a ID channel in heterostructures, a "quantum 
wire" || , the interactions are usually screened due to the 
presence of metallic gates near the channel. This leads to 
a Luttinger liquid || characterized by an interaction con- 
stant g (we only discuss the spinless single-channel case 
in the following). The noninteracting case corresponds 
to g = I, and the presence of (repulsive) Coulomb inter- 
actions implies g < 1. Quantities of principal interest are 
the channel conductance in the presence of interactions 
and the capacitance C = Q/U of a junction or impurity 
(U is the two-terminal voltage, and Q denotes the charge 
on the junction). We present an approach that allows one 
to address these problems for finite voltage at arbitrary 
interaction strength and junction transmittance. 

Our treatment is based on the standard bosonization 
approach jl] ||, which is applicable in the low-energy 
regime where only excitations near the Fermi surface are 
relevant. The electron creation operator can be expressed 
in terms of boson phase fields 9(x) and 4>(x), 



exp [ipkFX + iy/Tr[p9(x) + 4>( x )}] 

P =± 



where lj c — vpkp is the proper bandwidth cutoff for the 
linearized dispersion relation employed in the bosoniza- 
tion (we put h = 1). The phase fields obey the equal-time 
commutation relations 



[<P(x),6(x')] = -(i/2)sgn(x - x') , 

such that the canonical momentum for the 9 field is II = 
d x 4>- The boson representation for the electron density 
operator is then given by 



p{x) 



fej? 

7T 



-.d x 9(x) + — cos[2fc F ir + 2^9{x)] 

/IT 7T 



(1) 



The first term is the background charge, the second term 
stands for the sum of right- and left-moving densities p± , 
and the last term describes interference between right- 
and left-movers ||]. The clean Luttinger liquid is de- 
scribed by the Euclidean action M 



Sr, = — ^- / dxdr 



±(d T 9) 2 + ±(d x 9) 2 



(2) 



where vf is the Fermi velocity and g < 1 the interaction 
constant. A short-ranged impurity at x = results in 
the generic contribution @] 



S i = v J dTCOs[2^6(0,T)\ 



(3) 



where the dimensionless impurity strength A = irV/u c 
tunes the junction resistance. 

Now let us consider a ID quantum wire coupled to 
external reservoirs, see Fig. 1. The coupling of the 
ID channel to the 2D or 3D reservoirs is assumed 
to occur by adiabatic widening of the channel. This 
models ideal reservoirs in the sense of Landauer [Q. 
The reservoirs are macroscopic, totally incoherent de- 
vices kept at fixed chemical potential, and there are 
no reflections of particles entering the reservoirs. The 
left (right) reservoir at i — * — oo (x —> oo) has 
chemical potential /i_oo (poo), and the difference in 
chemical potentials is the two-terminal voltage U = 
(fi-oo — Poo)/e- The mean chemical potential (p-oo + 
p oa )/2 corresponds to the background charge kpj'K. 
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FIG. 1. General setup studied in this work: A ID quan- 
tum wire is coupled to 3D reservoirs, which are held at 
chemical potentials p_oc and ^oo- The left reservoir injects 
right-movers at the stationary excess density (p+), and the 
right one injects left-movers with density The striped 

area stands for a scatterer at x = 0. 

Applying a voltage means that one has a non-equilibrium 
influx of currents from the reservoirs: The reservoir at 
x — > — oo injects right-movers at some stationary excess 
density and similarly the density of left-movers at x — > oo 
is diminished, 



(p±(x))=± 



eU 



as x 



=Foo 



(4) 



While in a clean system the densities of right- and left- 
movers imposed by the reservoirs spread homogeneously 
along the wire, in the presence of scatterers one does 
not know beforehand what the densities of right- (left-) 
movers at x — > oo (x — » — oo) are. These densities follow 
from our theory. 

The boundary condition formulated in terms of 
the average densities of right- or left-moving fcrmions in 
the quantum wire as the reservoirs are approached, can 
equivalently be expressed in terms of the 9 field by noting 
that 



P+ + P- 



-.d x e 



The boundary condition for the boson phase field is there- 
fore 



OX Vf ot 



Too,r = 0)) = 



eU 



2^/ttvf 



(5) 



where the time r = has been picked by convention. 

To study the inhomogeneous quantum wire in presence 
of an external voltage, we consider the generating func- 
tional Z(y,fi) — (ex-p[2y/nifi9(y)]). We formally solve 
for Z by introducing an auxiliary field q(r) = 2^/n 9(x = 
0, r), with the constraint enforced by a Lagrange multi- 
plier field A(t). Then one has the effective action 



S e [6,A,q] 



VF 
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dxdr 



-{d T ef 



+ V J dr cosq(r) - 2 % Ari / u%,0) (6) 
+ i I dr A(t) [2^(9(0, r) - q(r)] . 



The 9 part of this effective action is Gaussian and can 
therefore be treated exactly by solving the classical Euler- 
Lagrange equations. One can always decompose the so- 
lution 6 into an homogeneous part Oh for the equilibrium 
system (U = 0), and a particular solution 9 p subject to 
the boundary condition (^). At the same time, we de- 
compose the field A — Ah 



A p , such that 



1 d 2 9 p 1 d 2 9 p 2y/^i 



v F dr 2 



g 2 dx 2 



vf 



S(x)A p (t) 



(J) 



The most general solution permitted by Eq. (R) which 
fulfills Eq. (H) requires a r independent A p and takes the 
form 

^r) = ^-^-\x\-ir e -^. (8) 

2y/TT 2y/TTVF 2^/TT 

The quantity tp is related to the zero mode of the La- 
grange multiplier field, A p = ieip/2irg 2 . 
With the boson propagators 



F(x,u) 



exp(-\gwx/v F \) 



the homogeneous part can be written in terms of the 
Fourier components A/ l (w), 



du> 
2^ 



[A h (u)F(x, u) - pF(x - y, u)]. 



Inserting 9^ + 9 p into Eq. (|6|), the action becomes Gaus- 
sian in Ah, which can therefore easily be integrated out. 
Adding the particular solution (0) also onto g(r) in view 
of q(r) = 2y / 7r6>(0, t), one obtains for the generating func- 
tional 



(9) 



(exp[2V^ip9(x)}) = W{xY / e -</*«¥>M/« 

f du ,F(x,uj) 
/ 2^ M 



x exp 



ifi q 



F(0,u>) 



where the average over the zero modes qo and <p of 
the auxiliary fields and over the q fluctuations q(r) = 
(27r) _1 J dujq(ui) exp(zwT) has to be taken using the ac- 
tion 



S 



cLo q(u))q(—u>) eip 



dr q(r) 



2vr 4F(0,w) 2ng 2 
+ V J dr cos[<jo - ie(U - ip)r + q(r)] . (10) 

The function W(x) — (1 + \x\/a>)~ 9 with the microscopic 
lengthscale a — vp/2guj c does not depend on impurity 
properties. 
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The effect of the external voltage can now be read off 
from Eqs. (g) and @. First, the average density [the 
first two terms in Eq. ([!])] is in general discontinuous at 
the impurity location due to particle reflection, 



9 
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sgn(x) 



(11) 



For the noninteracting case, g = 1, one can show from 
the exact solution of the equivalent Schrodinger equation 
that ip is the usual four-terminal voltage measured near 
the barrier ||. In that case, ip = RU, where R = 1/(1 + 
A~ 2 ) is the reflection coefficient of the barrier. 

Second, in the absence of a scatterer, A = 0, one always 
finds if = (see below), and the right- and left-moving 
densities are spatially homogeneous along the wire. Since 
they are determined by Eq. (0), the current is 



d9\ 
dr) 



ev F ({p+) - (p_)) 



U 



This yields the perfect two-terminal conductance e 2 /h in 
agreement with recent theoretical work [ [Tof and an exper- 
imental study of a quantum wire |7j. The conductance 
ge 2 /h discussed in Ref. fa] is not the two-terminal con- 
ductance but a low-frequency microwave conductance. 

To describe coupling to an external voltage, previous 
studies have often added a term to the Hamiltonian of 
the form |Q 



(12) 



where (p is the "voltage drop" |12|] . From Eq. (|10|) , if one 
tentatively identifies ip with (p, one observes that Eq. Jl2] ) 
should be modified by a factor l/g 2 . This factor can be 
understood in terms of the interaction energy of density 
fluctuations with the nonequilibrium background charge 
(O) deviating from kp/n, 



,/.,■ I dx' -^=d x 6{x) UJx - x') ( e</7Sgn(x ' ) 
Jtt V 2ixvp 



20? W 



dx sgn(x)d x 9(x) , 



where U c (x— x') is the screened Coulomb interaction, and 
the last line holds for a Luttinger liquid. Furthermore, it 
should be noted that in Eq. ( jl(i| ) the cosine part due to 
the impurity has acquired a term linear in time, which 
is reminiscent of the Josephson relations. This shows 
that the external voltage cannot be fully incorporated by 
simply adding a term like Eq. (|l^) to the Hamiltonian. 
In general, it is necessary to treat external reservoirs via 
boundary conditions ]Tq| . 

From the above considerations, we see that the action 
( |To| ) describes a voltage-biased ID quantum wire contain- 
ing a scatterer for the entire range of parameters. When 
evaluating Eq. (0), we still have to average over the zero 



mode of the Lagrange multiplier field. Hence ip is gener- 
ally a fluctuating quantity. In the two limiting cases of 
transmission zero and one, the fluctuations do vanish. In 
the latter case, A = 0, the q average is Gaussian, and one 
finds <p = due to the infrared singularity of the first 
term in the action @]. On the other hand, for A — ► oo, 
the cosine term in Eq. ( [To| ) strictly enforces ip = U. In 
the following, we will only consider the two fixed-point 
values ip — and U corresponding to a very small and a 
very high barrier, respectively. Near these limiting cases, 
further analytic progress can be made. 

We start by calculating the nonequilibrium electron 
density. Equivalent to an explicit real-time calculation, 
we first analytically continue to imaginary values of U 
and ip, and after performing the q average, we rotate back 
to real values of U and p. Let us first study the case of 
a very weak scatterer, such that we can put tp = and 
then use perturbation theory in the impurity strength. 
The antisymmetric charge distribution 

q(x) = -e[(p(x)) - (p(-x))}/2 (x > 0) 

can be computed by expressing the density operator (Q) 
in terms of the generating functional (g). Lowest-order 
perturbation theory in A yields 

q[x) = sin(2fcFx) ^L (a ./ a) -Gr-Va> (13) 

x [eU/oJcf- 1 ' 2 J g _ 1/2 (geUx/v F ) , 

where J u (x) is a Bessel function of the first kind [ pT[ . 
Near the impurity, for x <C (ell /vf) -1 , from properties 
of the Bessel function, 



q{x) 



eXkp 



2^r(g)T(g + 1/2) 



sm(2kpx) (eU/2u> ( 



The asymmetric charge mode q(x) is 2fcp-periodic but 
decays only on the lengthscale vp/eU. There is no lo- 
calized charge sitting on the impurity. However, the re- 
integration over ( p"3| ) gives a finite total charge, and hence 
a finite nonlinear capacitance C = Q/U, 



C(U) 



e 2 \/uj c 



eU 



V^r( 3 )r( 5 + i/2) V2^. 



(14) 



This lowest-order perturbational result in the weak- 
scattering regime breaks down for small voltages, 
et//2w c « AV(2-2s). 

In the opposite case of a strong scatterer, A > 1, we can 
put tp = U . From Eq. (^|), the 2k F -paxt of the asymmetric 
charge mode takes the form 

Q2k F (x) — — cos(2kpx) s\Yi(eUx/vp){x/ o)~ 9 , 

which implies a finite total charge Q2k F - This charge 
turns out to be linear in U , and therefore one has a finite 
U = capacitance 
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2k r 



-g-°r(2-g)sm -(1 - g) 



(15) 



However, qih F (x) does not include the constant term q = 
el) jlixvF coming from the change in background density 
(pd|). That term leads to a charge Qo = Le 2 U '/ 'Attvf, 
which diverges with the system length L. The applied 
voltage polarizes the capacitance between the wire and 
the metallic screening gate. This large shunt capacitance 
renders the observation of charging effects, i.e. of C2k F , 
impossible for a single scatterer. A similar situation is en- 
countered for a single tunnel junction coupled to metallic 
leads, where charging effects are normally absent 15 1. 

For an island formed by two strong impurities, the 2kp 
capacitance is observable since the capacitance of the is- 
land is not affected by the shunt capacitance between the 
remaining wire and the gate. Taking two impurities at 
x = ±R/2 and applying our boundary condition (||), we 
now have two zero modes from the respective Lagrange 
multipliers. One (tp) corresponds to the four-terminal 
voltage found in the single-impurity case, and the other 
{f G ) corresponds to a gate voltage applied to the island 
—R/2 < x < R/2. While again analytic results are not 
available for the entire range of parameters, it is possible 
to calculate the total charge Qi sitting on the island in 
the limit of large barriers and for UfR ^> 1- We obtain 



Qi 



eR 



(k F + eip G /2v F ) + 2C 2 k F <P G 



er(l- ff )sin[7r(l- 5 )/2]/7nf 



where C2k F is given in Eq. (|15|). The first term ~ R 
arises due to the slow component of the density operator 
(|l|), while the remaining two terms come from the 2k F 
component. The resonant tunneling condition can then 
be derived by noting that Qi/e is confined to integer 
values in the large-barrier limit under consideration here. 
Hence, the spacing of the resonances as a function of ip G 
is found to be 



eAip G = (R/2nv F + 2C 2kF /e 2 y 



(17) 



Therefore C2k F leads to an experimentally measurable 
resonance shift compared to previous results j5| which 
neglected charging effects in interacting ID metals. 

Let us finally comment on the issue of electroneutrality 
in a Luttinger liquid. The spatial change in the back- 
ground charge density ([ll]) induces an influence charge 
density of opposite sign on the metallic gate, such that 
overall electroneutrality is maintained. However, the 
Luttinger liquid interaction is not able to enforce elec- 
troneutrality within the ID quantum wire alone. In 
the complete absence of a gate, since there is no in- 
ternal screening within the wire, the Coulomb potential 
becomes long-ranged, U c (x — x') ~ \x — x'| _1 . From 
the Euler-Lagrange equation corresponding to Eq. (0) 
and the boundary condition (0), one finds immediately 



that then ip — 0. As a consequence, for a long-ranged 
Coulomb potential, electroneutrality is maintained auto- 
matically within the wire. 

In conclusion, the theoretical description of a voltage- 
biased ID quantum wire (Luttinger liquid) containing 
elastic potential scatterers has been given. Our boundary 
condition method can easily be adopted to the case of 
more than one channel (e.g., the spin-^ case), and it can 
straightforwardly be generalized to a real-time and finite- 
temperature formalism. 
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